Rockfish Historical Harvest Estimates

This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.

Issue	Howard	Bayes
Time series	1999 - present	1977 - present
Bias in SWHS	Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests.	Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases	Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data.	Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations	Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas.	Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation	Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with.	By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Results

• Model: LB_fit_3pH

• Model Run: *Gen4int_indcomp_swhsR_FULL_pHB4pars_re0full_altwt_slepsdwn_thru2024_3e+05__2025-11-06* (Naming convention: modelname_thruYEAR_iters_rundate)

• Proportion of parmaters converged (Rhat <= 1.1): 0.9958

• See end of document for summary of unconverged parameters.

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

And black rockfish releases…

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Total Biomass Removal Estimates

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Model fit

Logbook residuals

Figure 20.- Residuals from logbook harvests.

SWHS residuals

Figure 21.- Residuals from SWHS harvests.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

Figure 23.- Mean percent of harvest by charter anglers.

When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(slope|non-pelagic & non-yellowye) For release estimates

## NULL

Weight Fits

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.5204082

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.1 while 0.9 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable	area_n	year	user	mean	median	lower	upper	Rhat	sd	skew	area
p_pelagic	13	1	1	0.301	0.277	0.103	0.616	1.144	0.135	0.178	NSEI
p_pelagic	13	1	2	0.124	0.105	0.033	0.318	1.131	0.075	0.255	NSEI
p_pelagic	13	10	1	0.317	0.299	0.110	0.615	1.118	0.133	0.136	NSEI
p_pelagic	13	10	2	0.132	0.115	0.036	0.330	1.138	0.076	0.225	NSEI
p_pelagic	13	11	1	0.314	0.295	0.113	0.605	1.113	0.132	0.144	NSEI
p_pelagic	13	11	2	0.135	0.114	0.038	0.333	1.147	0.078	0.268	NSEI
p_pelagic	13	12	1	0.326	0.310	0.117	0.620	1.108	0.134	0.119	NSEI
p_pelagic	13	12	2	0.135	0.117	0.039	0.339	1.163	0.078	0.232	NSEI
p_pelagic	13	13	1	0.326	0.311	0.118	0.607	1.113	0.131	0.115	NSEI
p_pelagic	13	13	2	0.137	0.120	0.039	0.330	1.137	0.076	0.225	NSEI
p_pelagic	13	14	1	0.309	0.292	0.121	0.587	1.098	0.121	0.140	NSEI
p_pelagic	13	14	2	0.135	0.120	0.040	0.314	1.124	0.073	0.206	NSEI
p_pelagic	13	15	1	0.283	0.267	0.102	0.557	1.127	0.120	0.134	NSEI
p_pelagic	13	15	2	0.130	0.115	0.040	0.312	1.124	0.070	0.214	NSEI
p_pelagic	13	16	1	0.321	0.302	0.128	0.597	1.109	0.124	0.153	NSEI
p_pelagic	13	16	2	0.143	0.128	0.044	0.325	1.117	0.074	0.203	NSEI
p_pelagic	13	17	1	0.334	0.321	0.139	0.595	1.080	0.120	0.108	NSEI
p_pelagic	13	17	2	0.149	0.133	0.047	0.334	1.088	0.074	0.216	NSEI
p_pelagic	13	18	1	0.345	0.331	0.144	0.611	1.073	0.121	0.116	NSEI
p_pelagic	13	18	2	0.155	0.141	0.050	0.343	1.064	0.076	0.184	NSEI
p_pelagic	13	19	1	0.347	0.332	0.144	0.625	1.051	0.124	0.121	NSEI
p_pelagic	13	19	2	0.155	0.142	0.053	0.336	1.057	0.074	0.176	NSEI
p_pelagic	13	2	1	0.300	0.278	0.099	0.603	1.126	0.135	0.163	NSEI
p_pelagic	13	2	2	0.126	0.106	0.035	0.314	1.135	0.074	0.269	NSEI
p_pelagic	13	20	1	0.357	0.345	0.163	0.616	1.058	0.119	0.101	NSEI
p_pelagic	13	20	2	0.164	0.150	0.054	0.355	1.058	0.077	0.182	NSEI
p_pelagic	13	21	1	0.374	0.364	0.167	0.645	1.045	0.125	0.080	NSEI
p_pelagic	13	21	2	0.173	0.160	0.059	0.371	1.039	0.080	0.162	NSEI
p_pelagic	13	22	2	0.172	0.159	0.065	0.354	1.026	0.075	0.173	NSEI
p_pelagic	13	23	2	0.111	0.101	0.041	0.237	1.017	0.051	0.197	NSEI
p_pelagic	13	3	1	0.298	0.275	0.097	0.606	1.161	0.137	0.168	NSEI
p_pelagic	13	3	2	0.125	0.106	0.034	0.326	1.189	0.077	0.247	NSEI
p_pelagic	13	4	1	0.306	0.286	0.101	0.619	1.151	0.135	0.148	NSEI
p_pelagic	13	4	2	0.124	0.103	0.034	0.331	1.153	0.077	0.274	NSEI
p_pelagic	13	5	1	0.301	0.279	0.097	0.607	1.151	0.134	0.164	NSEI
p_pelagic	13	5	2	0.127	0.109	0.034	0.321	1.163	0.076	0.236	NSEI
p_pelagic	13	6	1	0.307	0.288	0.102	0.615	1.146	0.135	0.141	NSEI
p_pelagic	13	6	2	0.126	0.108	0.033	0.320	1.180	0.075	0.239	NSEI
p_pelagic	13	7	1	0.307	0.288	0.109	0.607	1.137	0.133	0.143	NSEI
p_pelagic	13	7	2	0.128	0.108	0.034	0.326	1.140	0.076	0.264	NSEI
p_pelagic	13	8	1	0.308	0.291	0.107	0.601	1.127	0.132	0.129	NSEI
p_pelagic	13	8	2	0.128	0.111	0.034	0.315	1.152	0.073	0.232	NSEI
p_pelagic	13	9	1	0.309	0.290	0.107	0.605	1.124	0.131	0.145	NSEI
p_pelagic	13	9	2	0.128	0.110	0.034	0.334	1.158	0.077	0.235	NSEI
p_pelagic	14	2	2	0.410	0.404	0.184	0.667	1.011	0.125	0.048	NSEO
p_black	14	3	1	0.918	0.945	0.688	0.993	1.017	0.083	-0.327	NSEO
p_pelagic	16	1	1	0.345	0.332	0.116	0.642	1.225	0.140	0.093	SSEO
p_pelagic	16	1	2	0.326	0.311	0.115	0.625	1.219	0.136	0.110	SSEO
p_pelagic	16	10	1	0.350	0.342	0.122	0.653	1.194	0.139	0.057	SSEO
p_pelagic	16	10	2	0.333	0.320	0.118	0.624	1.213	0.136	0.096	SSEO
p_pelagic	16	11	1	0.353	0.341	0.128	0.646	1.160	0.139	0.086	SSEO
p_pelagic	16	11	2	0.332	0.320	0.118	0.605	1.203	0.131	0.091	SSEO
p_pelagic	16	12	1	0.355	0.345	0.128	0.650	1.177	0.138	0.072	SSEO
p_pelagic	16	12	2	0.331	0.315	0.121	0.614	1.163	0.134	0.120	SSEO
p_pelagic	16	13	1	0.357	0.343	0.132	0.640	1.162	0.136	0.103	SSEO
p_pelagic	16	13	2	0.338	0.324	0.124	0.623	1.165	0.133	0.105	SSEO
p_pelagic	16	14	1	0.345	0.332	0.132	0.621	1.144	0.131	0.099	SSEO
p_pelagic	16	14	2	0.337	0.323	0.123	0.625	1.135	0.132	0.106	SSEO
p_pelagic	16	15	1	0.351	0.339	0.139	0.623	1.153	0.130	0.092	SSEO
p_pelagic	16	15	2	0.340	0.327	0.132	0.619	1.127	0.129	0.101	SSEO
p_pelagic	16	16	1	0.349	0.336	0.136	0.627	1.133	0.131	0.099	SSEO
p_pelagic	16	16	2	0.342	0.332	0.134	0.609	1.122	0.127	0.079	SSEO
p_pelagic	16	17	1	0.404	0.394	0.181	0.669	1.095	0.129	0.078	SSEO
p_pelagic	16	17	2	0.360	0.348	0.135	0.641	1.127	0.133	0.090	SSEO
p_pelagic	16	18	1	0.390	0.384	0.163	0.660	1.096	0.130	0.046	SSEO
p_pelagic	16	18	2	0.360	0.349	0.152	0.629	1.106	0.128	0.086	SSEO
p_pelagic	16	19	1	0.379	0.367	0.163	0.651	1.093	0.128	0.093	SSEO
p_pelagic	16	19	2	0.363	0.354	0.145	0.632	1.077	0.130	0.069	SSEO
p_pelagic	16	2	1	0.342	0.326	0.122	0.645	1.237	0.141	0.113	SSEO
p_pelagic	16	2	2	0.324	0.307	0.104	0.625	1.253	0.137	0.124	SSEO
p_pelagic	16	20	1	0.387	0.378	0.172	0.650	1.085	0.125	0.072	SSEO
p_pelagic	16	20	2	0.369	0.361	0.153	0.628	1.076	0.125	0.064	SSEO
p_pelagic	16	21	1	0.371	0.364	0.162	0.622	1.057	0.122	0.058	SSEO
p_pelagic	16	21	2	0.369	0.357	0.158	0.631	1.047	0.125	0.096	SSEO
p_pelagic	16	22	2	0.377	0.371	0.166	0.638	1.035	0.123	0.049	SSEO
p_pelagic	16	23	2	0.300	0.287	0.119	0.547	1.023	0.112	0.116	SSEO
p_pelagic	16	3	1	0.348	0.327	0.124	0.646	1.246	0.142	0.148	SSEO
p_pelagic	16	3	2	0.329	0.316	0.107	0.634	1.240	0.140	0.093	SSEO
p_dsr	16	30	1	0.419	0.419	0.393	0.447	1.015	0.014	0.000	SSEO
p_pelagic	16	4	1	0.344	0.329	0.122	0.632	1.245	0.140	0.107	SSEO
p_pelagic	16	4	2	0.325	0.311	0.105	0.626	1.202	0.136	0.103	SSEO
p_pelagic	16	5	1	0.347	0.333	0.123	0.639	1.214	0.140	0.100	SSEO
p_pelagic	16	5	2	0.325	0.308	0.107	0.611	1.190	0.136	0.125	SSEO
p_pelagic	16	6	1	0.346	0.330	0.118	0.634	1.208	0.141	0.114	SSEO
p_pelagic	16	6	2	0.326	0.317	0.100	0.615	1.185	0.136	0.066	SSEO
p_pelagic	16	7	1	0.345	0.332	0.119	0.631	1.204	0.138	0.094	SSEO
p_pelagic	16	7	2	0.331	0.314	0.117	0.619	1.223	0.135	0.126	SSEO
p_pelagic	16	8	1	0.347	0.332	0.124	0.639	1.214	0.138	0.109	SSEO
p_pelagic	16	8	2	0.330	0.315	0.110	0.626	1.195	0.137	0.110	SSEO
p_pelagic	16	9	1	0.352	0.341	0.120	0.649	1.187	0.140	0.079	SSEO
p_pelagic	16	9	2	0.329	0.313	0.116	0.618	1.193	0.133	0.121	SSEO
p_pelagic	3	1	1	0.546	0.558	0.262	0.780	1.011	0.135	-0.089	PWSI
p_black	8	40	1	0.523	0.422	0.341	0.969	1.143	0.220	0.459	afognak
p_black	8	48	1	0.537	0.493	0.411	0.969	1.043	0.148	0.298	afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable	area_n	year/species	mean	median	lower	upper	Rhat	sd	skew	year	area
beta2_pH	10	1	1.579	0.895	0.320	7.822	1.145	1.878	0.364		northeast
beta0_pH	11	4	0.860	0.897	-0.343	1.858	1.109	0.537	-0.069		CSEO
beta1_pH	11	1	3.488	3.441	2.436	4.653	1.123	0.568	0.083		CSEO
beta1_pH	11	4	1.926	1.918	0.786	3.104	1.103	0.579	0.014		CSEO
beta1_pelagic	11		1.185	1.158	0.711	1.824	1.142	0.263	0.103		CSEO
beta2_yellow	11		-6.842	-5.681	-20.016	-0.411	1.118	5.237	-0.222		CSEO
beta0_pH	13	4	0.478	0.588	-1.061	1.576	1.179	0.688	-0.160		NSEI
beta1_pH	13	4	1.494	1.401	0.346	3.045	1.178	0.688	0.135		NSEI
beta0_pelagic	13		-1.008	-1.103	-1.713	-0.004	1.341	0.475	0.200		NSEI
beta1_pelagic	13		2.810	2.891	1.477	3.843	1.325	0.614	-0.132		NSEI
beta2_pelagic	13		0.188	0.162	0.098	0.443	1.263	0.092	0.283		NSEI
beta3_pelagic	13		31.164	30.680	24.723	38.969	1.348	3.996	0.121		NSEI
beta0_pH	15	4	0.667	0.684	-0.291	1.576	1.561	0.459	-0.037		SSEI
beta1_pH	15	4	1.345	1.343	0.290	2.334	1.514	0.521	0.004		SSEI
beta1_pelagic	15		0.090	0.013	0.000	0.603	1.108	0.162	0.477		SSEI
beta0_pelagic	16		-0.719	-0.733	-1.524	0.026	1.695	0.455	0.031		SSEO
beta1_pelagic	16		1.858	1.878	0.827	3.034	1.663	0.629	-0.032		SSEO
beta2_pelagic	16		0.560	0.271	0.108	3.745	1.195	0.875	0.330		SSEO
beta3_pelagic	16		29.877	30.135	24.912	34.665	1.184	2.591	-0.100		SSEO
mu_beta2_pH	2	1	1.583	0.882	0.198	8.069	1.161	2.029	0.346		NG
beta0_pH	2	1	1.418	1.426	1.011	1.723	1.121	0.180	-0.044		NG
beta1_pH	2	1	1.545	1.535	1.037	2.134	1.104	0.282	0.036		NG
mu_beta0_pH	3	4	0.688	0.729	-0.392	1.582	1.141	0.512	-0.080		PWSI
tau_beta1_pH	3	4	5.104	3.020	0.247	21.889	1.144	7.293	0.286		PWSI
tau_beta2_pH	3	1	4.714	1.092	0.044	39.204	1.112	10.952	0.331		PWSI
beta0_pH	3	1	1.354	1.403	0.366	1.768	1.236	0.303	-0.162		PWSI
beta1_pH	3	1	1.865	1.785	1.108	3.437	1.207	0.514	0.156		PWSI
beta2_pH	6	1	1.361	0.523	0.126	8.858	1.146	2.224	0.377		SOKO2SAP
beta0_pH	8	1	-0.560	-0.539	-1.180	-0.033	1.175	0.294	-0.071		afognak
beta2_pH	8	1	1.024	0.402	0.150	6.871	1.256	1.714	0.363		afognak
beta3_pH	8	1	35.023	34.860	29.093	42.115	1.251	3.614	0.045		afognak
beta2_pH	9	1	1.365	0.610	0.158	8.275	1.112	1.990	0.379		eastside
lambda_H	15		0.019	0.014	0.003	0.057	1.249	0.047	0.107		SSEI
lambda_H	6		3.837	0.117	0.008	32.088	1.144	9.178	0.405		SOKO2SAP

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.01	636.00	292.00	1,186.00	358.00
Prop. of total	0.69	0.31	0.77	0.23
Prop. median less than 1	0.21	0.16	0.31	0.41
Prop. median less than 10	0.33	0.41	0.48	0.70
Prop. median less than 50	0.49	0.62	0.65	0.77
Median of the median values	53.28	21.26	11.16	1.88
Mean of the median values	584.32	627.73	332.50	261.96
Mean skew of posteriors	0.23	0.29	0.27	0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.1	31.00	3.00	68.00	35.00
Prop. of total	0.91	0.09	0.66	0.34
Prop. median less than 1	0.13	1.00	0.71	0.63
Prop. median less than 10	0.16	1.00	0.76	0.91
Prop. median less than 50	0.19	1.00	0.82	0.91
Median of the median values	295.80	0.09	0.01	0.28
Mean of the median values	443.61	0.36	51.42	14.85
Mean skew of posteriors	0.25	0.21	0.16	0.16

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable	Freq
beta0_pelagic	2
beta0_pH	5
beta1_pelagic	3
beta1_pH	4
beta2_pelagic	2
beta2_pH	4
beta2_yellow	1
beta3_pelagic	2
beta3_pH	1
Bp_ay	5
Bp_ayu	4
eps_bl	1
eps_pel	6
eps_pH	11
H_ay	2
H_ayg	1
H_ayu	1
Hb_ay	6
Hb_ayg	2
Hb_ayu	3
Ho_ay	2
Ho_ayg	2
Ho_ayu	3
Hp_ay	12
Hp_ayg	3
Hp_ayu	7
Hs_ayu	1
Htrend_ay	4
Hy_ayu	1
lambda_H	2
mu_beta0_pH	1
mu_beta2_pH	1
p_black	1
p_pelagic	62
p_slope	1
pH	5
R_ay	2
R_ayg	6
R_ayu	8
Rb_ay	1
Rb_ayg	5
Rb_ayu	1
Rd_ay	1
Rd_ayg	2
Rd_ayu	1
re_black	1
re_pelagic	19
re_pH	16
Ro_ay	18
Ro_ayg	13
Ro_ayu	33
Rp_ay	2
Rp_ayg	6
Rp_ayu	2
Rs_ay	3
Rs_ayg	9
Rs_ayu	16
Ry_ay	11
Ry_ayg	4
Ry_ayu	23
tau_beta1_pH	1
tau_beta2_pH	1

Parameter estimates:

Summary Table of Parameter Estimates

Parameter	mean	sd	Lower_CI	Median	Upper_CI
mu_bc_H[1]	-0.020	0.067	-0.145	-0.023	0.124
mu_bc_H[2]	0.078	0.054	-0.018	0.076	0.190
mu_bc_H[3]	-0.267	0.068	-0.395	-0.266	-0.134
mu_bc_H[4]	-0.921	0.179	-1.280	-0.921	-0.576
mu_bc_H[5]	1.308	1.064	0.140	1.071	3.942
mu_bc_H[6]	-2.376	0.313	-2.999	-2.372	-1.773
mu_bc_H[7]	-0.456	0.103	-0.665	-0.455	-0.262
mu_bc_H[8]	0.253	0.356	-0.326	0.222	0.984
mu_bc_H[9]	-0.244	0.126	-0.481	-0.248	0.015
mu_bc_H[10]	0.058	0.067	-0.069	0.056	0.198
mu_bc_H[11]	-0.124	0.037	-0.197	-0.124	-0.050
mu_bc_H[12]	-0.254	0.105	-0.476	-0.251	-0.053
mu_bc_H[13]	-0.157	0.072	-0.290	-0.159	-0.020
mu_bc_H[14]	-0.311	0.087	-0.491	-0.309	-0.147
mu_bc_H[15]	-0.338	0.048	-0.433	-0.336	-0.242
mu_bc_H[16]	-0.161	0.371	-0.805	-0.186	0.649
mu_bc_R[1]	1.041	0.125	0.799	1.040	1.279
mu_bc_R[2]	1.014	0.076	0.854	1.017	1.156
mu_bc_R[3]	0.226	0.092	0.044	0.226	0.404
mu_bc_R[4]	0.676	0.197	0.271	0.683	1.049
mu_bc_R[5]	1.686	0.367	0.970	1.680	2.433
mu_bc_R[6]	-1.404	0.349	-2.096	-1.413	-0.722
mu_bc_R[7]	0.505	0.169	0.158	0.509	0.826
mu_bc_R[8]	0.657	0.149	0.378	0.657	0.954
mu_bc_R[9]	0.403	0.163	0.062	0.409	0.708
mu_bc_R[10]	0.899	0.105	0.693	0.901	1.105
mu_bc_R[11]	1.106	0.074	0.962	1.105	1.257
mu_bc_R[12]	0.874	0.165	0.564	0.872	1.196
mu_bc_R[13]	1.037	0.091	0.860	1.038	1.212
mu_bc_R[14]	0.855	0.119	0.624	0.857	1.090
mu_bc_R[15]	0.798	0.095	0.606	0.799	0.988
mu_bc_R[16]	1.357	0.115	1.126	1.361	1.582
tau_pH[1]	2.432	0.182	2.090	2.423	2.794
tau_pH[2]	1.611	0.150	1.337	1.610	1.914
tau_pH[3]	1.194	0.110	0.993	1.191	1.416
tau_pH[4]	1.015	0.138	0.770	1.007	1.307
tau_pH[5]	0.427	0.076	0.297	0.420	0.592
beta0_pH[1,1]	0.581	0.171	0.219	0.587	0.898
beta0_pH[2,1]	1.418	0.180	1.011	1.426	1.723
beta0_pH[3,1]	1.354	0.303	0.366	1.403	1.768
beta0_pH[4,1]	1.572	0.251	0.997	1.593	2.004
beta0_pH[5,1]	-0.735	0.376	-1.589	-0.710	-0.092
beta0_pH[6,1]	-0.476	0.452	-1.732	-0.354	0.090
beta0_pH[7,1]	0.757	0.011	0.736	0.757	0.778
beta0_pH[8,1]	-0.560	0.294	-1.180	-0.539	-0.033
beta0_pH[9,1]	-0.249	0.279	-0.868	-0.212	0.190
beta0_pH[10,1]	0.278	0.190	-0.094	0.287	0.653
beta0_pH[11,1]	-0.032	0.337	-0.662	-0.062	0.719
beta0_pH[12,1]	0.632	0.198	0.216	0.642	1.000
beta0_pH[13,1]	-0.654	0.150	-0.955	-0.651	-0.376
beta0_pH[14,1]	0.004	0.176	-0.363	0.006	0.339
beta0_pH[15,1]	0.041	0.185	-0.348	0.051	0.374
beta0_pH[16,1]	-0.914	0.220	-1.346	-0.915	-0.469
beta0_pH[1,2]	2.563	0.084	2.404	2.560	2.734
beta0_pH[2,2]	3.079	0.022	3.036	3.080	3.123
beta0_pH[3,2]	3.626	0.042	3.544	3.625	3.711
beta0_pH[4,2]	2.937	0.019	2.901	2.936	2.974
beta0_pH[5,2]	4.528	1.488	2.713	4.229	8.362
beta0_pH[6,2]	3.320	0.163	3.005	3.316	3.642
beta0_pH[7,2]	2.033	0.049	1.941	2.032	2.133
beta0_pH[8,2]	2.820	0.065	2.697	2.819	2.953
beta0_pH[9,2]	3.954	0.210	3.569	3.947	4.382
beta0_pH[10,2]	3.996	0.130	3.746	3.993	4.258
beta0_pH[11,2]	1.885	0.009	1.869	1.885	1.903
beta0_pH[12,2]	1.796	0.069	1.662	1.794	1.936
beta0_pH[13,2]	2.496	0.026	2.446	2.496	2.548
beta0_pH[14,2]	1.803	0.024	1.759	1.803	1.850
beta0_pH[15,2]	2.545	0.017	2.512	2.545	2.578
beta0_pH[16,2]	2.765	0.024	2.718	2.765	2.812
beta0_pH[1,3]	1.075	0.401	0.151	1.146	1.679
beta0_pH[2,3]	2.126	0.018	2.092	2.126	2.160
beta0_pH[3,3]	2.692	0.024	2.645	2.692	2.740
beta0_pH[4,3]	2.980	0.034	2.915	2.980	3.046
beta0_pH[5,3]	3.495	2.442	0.911	2.854	9.990
beta0_pH[6,3]	1.382	0.381	0.721	1.362	2.187
beta0_pH[7,3]	0.859	0.051	0.760	0.861	0.958
beta0_pH[8,3]	0.249	0.301	-0.469	0.291	0.653
beta0_pH[9,3]	-2.299	1.062	-4.732	-2.185	-0.501
beta0_pH[10,3]	-0.427	0.835	-2.414	-0.327	0.891
beta0_pH[11,3]	1.157	0.486	0.132	1.175	2.049
beta0_pH[12,3]	1.009	0.651	-0.366	0.999	2.261
beta0_pH[13,3]	1.047	0.481	-0.004	1.072	1.909
beta0_pH[14,3]	0.620	0.413	-0.277	0.651	1.364
beta0_pH[15,3]	0.774	0.414	-0.087	0.790	1.528
beta0_pH[16,3]	1.185	0.634	-0.095	1.181	2.465
beta0_pH[11,4]	0.860	0.537	-0.343	0.897	1.858
beta0_pH[12,4]	0.883	0.716	-0.662	0.885	2.246
beta0_pH[13,4]	0.478	0.688	-1.061	0.588	1.576
beta0_pH[14,4]	0.392	0.466	-0.630	0.429	1.240
beta0_pH[15,4]	0.667	0.459	-0.291	0.684	1.576
beta0_pH[16,4]	0.855	0.632	-0.495	0.889	2.047
beta0_pH[11,5]	2.971	0.145	2.698	2.968	3.264
beta0_pH[12,5]	-0.920	0.518	-2.018	-0.897	0.014
beta0_pH[13,5]	4.902	0.422	4.144	4.885	5.794
beta0_pH[14,5]	3.894	0.364	3.233	3.881	4.669
beta0_pH[15,5]	2.285	0.124	2.049	2.282	2.540
beta0_pH[16,5]	0.253	0.065	0.161	0.240	0.410
beta1_pH[1,1]	2.476	0.464	1.704	2.420	3.557
beta1_pH[2,1]	1.545	0.282	1.037	1.535	2.134
beta1_pH[3,1]	1.865	0.514	1.108	1.785	3.437
beta1_pH[4,1]	1.893	0.417	1.214	1.838	2.882
beta1_pH[5,1]	2.242	0.484	1.407	2.202	3.346
beta1_pH[6,1]	2.027	0.757	1.065	1.823	4.104
beta1_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,1]	2.228	0.777	1.136	2.084	4.126
beta1_pH[9,1]	1.968	0.600	1.188	1.847	3.688
beta1_pH[10,1]	2.098	0.362	1.456	2.084	2.878
beta1_pH[11,1]	3.488	0.568	2.436	3.441	4.653
beta1_pH[12,1]	2.894	0.292	2.341	2.885	3.469
beta1_pH[13,1]	3.451	0.239	3.018	3.438	3.940
beta1_pH[14,1]	2.857	0.330	2.263	2.840	3.570
beta1_pH[15,1]	2.414	0.252	1.969	2.405	2.941
beta1_pH[16,1]	4.087	0.303	3.472	4.093	4.689
beta1_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[1,3]	1.969	1.174	0.049	1.835	5.116
beta1_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,3]	4.735	2.537	0.619	4.224	11.082
beta1_pH[9,3]	7.669	2.717	3.360	7.305	13.825
beta1_pH[10,3]	6.277	2.571	2.717	5.789	12.506
beta1_pH[11,3]	2.127	0.523	1.110	2.119	3.177
beta1_pH[12,3]	3.214	0.669	1.933	3.203	4.654
beta1_pH[13,3]	1.929	0.517	0.921	1.912	3.025
beta1_pH[14,3]	2.148	0.469	1.273	2.136	3.140
beta1_pH[15,3]	2.039	0.484	1.123	2.026	3.018
beta1_pH[16,3]	3.695	0.685	2.329	3.706	5.110
beta1_pH[11,4]	1.926	0.579	0.786	1.918	3.104
beta1_pH[12,4]	3.189	0.744	1.798	3.156	4.763
beta1_pH[13,4]	1.494	0.688	0.346	1.401	3.045
beta1_pH[14,4]	1.976	0.538	0.960	1.954	3.137
beta1_pH[15,4]	1.345	0.521	0.290	1.343	2.334
beta1_pH[16,4]	3.536	0.695	2.186	3.500	5.032
beta1_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,1]	1.379	0.890	0.344	1.085	3.377
beta2_pH[2,1]	1.767	0.929	0.386	1.619	3.430
beta2_pH[3,1]	1.383	0.863	0.159	1.274	2.946
beta2_pH[4,1]	1.171	0.797	0.221	0.932	2.895
beta2_pH[5,1]	2.298	2.452	0.293	1.549	9.821
beta2_pH[6,1]	1.361	2.224	0.126	0.523	8.858
beta2_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,1]	1.024	1.714	0.150	0.402	6.871
beta2_pH[9,1]	1.365	1.990	0.158	0.610	8.275
beta2_pH[10,1]	1.579	1.878	0.320	0.895	7.822
beta2_pH[11,1]	0.379	0.201	0.172	0.325	0.967
beta2_pH[12,1]	1.573	0.807	0.614	1.336	3.672
beta2_pH[13,1]	1.032	0.446	0.545	0.918	2.306
beta2_pH[14,1]	1.093	0.626	0.387	0.933	3.134
beta2_pH[15,1]	0.991	0.559	0.380	0.852	2.653
beta2_pH[16,1]	1.003	0.310	0.613	0.932	1.815
beta2_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,3]	6.654	6.578	0.154	4.905	22.856
beta2_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,3]	1.539	1.531	0.118	1.055	6.375
beta2_pH[9,3]	0.348	0.883	0.088	0.195	1.404
beta2_pH[10,3]	0.481	0.917	0.109	0.278	2.172
beta2_pH[11,3]	-6.286	4.856	-18.752	-4.947	-0.724
beta2_pH[12,3]	-6.934	4.737	-19.149	-5.890	-1.220
beta2_pH[13,3]	-6.427	4.722	-18.635	-5.015	-1.076
beta2_pH[14,3]	-6.393	4.628	-18.457	-4.979	-1.285
beta2_pH[15,3]	-6.533	4.667	-18.559	-5.273	-1.040
beta2_pH[16,3]	-6.522	4.577	-18.808	-5.156	-1.507
beta2_pH[11,4]	-5.273	3.948	-14.736	-4.077	-0.752
beta2_pH[12,4]	-5.660	4.015	-15.711	-4.542	-1.027
beta2_pH[13,4]	-5.380	3.967	-15.173	-4.143	-0.814
beta2_pH[14,4]	-5.521	3.910	-15.256	-4.414	-1.069
beta2_pH[15,4]	-5.570	3.974	-15.302	-4.470	-0.844
beta2_pH[16,4]	-5.610	3.748	-15.015	-4.512	-1.377
beta2_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta3_pH[1,1]	34.866	1.083	32.828	34.886	37.279
beta3_pH[2,1]	32.368	1.303	29.786	32.368	34.953
beta3_pH[3,1]	32.919	1.469	29.927	32.917	35.868
beta3_pH[4,1]	32.642	1.829	29.677	32.513	36.433
beta3_pH[5,1]	27.694	1.235	26.131	27.515	31.175
beta3_pH[6,1]	33.121	3.115	26.663	32.974	40.789
beta3_pH[7,1]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,1]	35.023	3.614	29.093	34.860	42.115
beta3_pH[9,1]	32.569	2.569	28.179	32.199	38.639
beta3_pH[10,1]	32.358	1.239	30.082	32.315	34.976
beta3_pH[11,1]	31.163	1.072	29.273	31.091	33.275
beta3_pH[12,1]	30.790	0.467	29.810	30.803	31.680
beta3_pH[13,1]	31.423	0.539	30.394	31.417	32.485
beta3_pH[14,1]	32.494	0.742	31.199	32.455	34.039
beta3_pH[15,1]	30.915	0.800	29.342	30.918	32.537
beta3_pH[16,1]	30.106	0.354	29.438	30.099	30.817
beta3_pH[1,2]	34.543	5.416	25.577	34.480	43.495
beta3_pH[2,2]	34.703	5.423	25.474	34.852	43.517
beta3_pH[3,2]	30.405	7.179	18.676	30.252	42.362
beta3_pH[4,2]	30.306	7.166	18.599	30.263	42.318
beta3_pH[5,2]	30.506	7.297	18.568	30.500	42.388
beta3_pH[6,2]	30.741	7.220	18.558	30.755	42.464
beta3_pH[7,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[9,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[10,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[11,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[12,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[13,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[14,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[15,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[16,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[1,3]	38.305	3.199	30.847	39.597	42.481
beta3_pH[2,3]	34.469	5.455	25.475	34.489	43.516
beta3_pH[3,3]	30.377	7.135	18.637	30.394	42.248
beta3_pH[4,3]	30.500	7.138	18.784	30.477	42.510
beta3_pH[5,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[6,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[7,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,3]	42.028	1.373	38.275	42.082	43.903
beta3_pH[9,3]	36.204	1.900	34.528	35.442	41.810
beta3_pH[10,3]	37.755	2.448	34.672	37.071	43.421
beta3_pH[11,3]	40.696	1.069	38.800	40.653	42.736
beta3_pH[12,3]	41.992	0.742	40.319	42.244	42.859
beta3_pH[13,3]	41.399	0.730	40.009	41.409	42.829
beta3_pH[14,3]	40.477	0.675	39.128	40.494	41.730
beta3_pH[15,3]	40.629	0.650	39.194	40.641	41.893
beta3_pH[16,3]	41.846	0.391	41.044	41.907	42.622
beta3_pH[11,4]	40.616	1.187	38.359	40.494	42.981
beta3_pH[12,4]	41.857	0.787	40.223	42.111	42.950
beta3_pH[13,4]	42.268	1.190	39.837	42.566	43.895
beta3_pH[14,4]	40.479	0.734	38.998	40.524	41.802
beta3_pH[15,4]	41.100	1.147	38.360	41.056	43.418
beta3_pH[16,4]	41.743	0.423	40.928	41.785	42.567
beta3_pH[11,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[12,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[13,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[14,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[15,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[16,5]	25.000	0.000	25.000	25.000	25.000
beta4_pH[1,1]	0.264	0.021	0.226	0.264	0.307
beta4_pH[2,1]	0.211	0.005	0.201	0.211	0.221
beta4_pH[3,1]	0.173	0.007	0.159	0.173	0.186
beta4_pH[4,1]	0.294	0.012	0.271	0.294	0.316
beta4_pH[5,1]	0.419	0.277	0.086	0.385	0.922
beta4_pH[6,1]	0.784	0.402	0.325	0.733	1.445
beta4_pH[7,1]	0.858	0.340	0.355	0.809	1.646
beta4_pH[8,1]	0.159	0.069	0.046	0.153	0.309
beta4_pH[9,1]	0.211	0.084	0.062	0.208	0.388
beta4_pH[10,1]	0.001	0.001	0.000	0.001	0.003
beta4_pH[11,1]	0.447	0.012	0.424	0.447	0.470
beta4_pH[12,1]	0.641	0.022	0.598	0.641	0.686
beta4_pH[13,1]	0.631	0.013	0.605	0.631	0.657
beta4_pH[14,1]	0.670	0.030	0.612	0.669	0.731
beta4_pH[15,1]	0.796	0.021	0.756	0.797	0.838
beta4_pH[16,1]	1.218	0.042	1.141	1.216	1.306
beta4_pH[1,2]	0.308	0.094	0.136	0.300	0.507
beta4_pH[2,2]	0.514	0.023	0.470	0.513	0.561
beta4_pH[3,2]	0.446	0.012	0.422	0.446	0.469
beta4_pH[4,2]	0.654	0.038	0.581	0.653	0.730
beta4_pH[5,2]	0.733	0.918	0.112	0.530	2.733
beta4_pH[6,2]	0.916	0.795	0.202	0.737	2.876
beta4_pH[7,2]	0.894	0.446	0.357	0.810	1.914
beta4_pH[8,2]	0.849	0.199	0.494	0.837	1.265
beta4_pH[9,2]	1.195	1.235	0.306	0.886	4.042
beta4_pH[10,2]	0.457	0.051	0.361	0.457	0.559
beta4_pH[11,2]	1.270	0.049	1.177	1.270	1.368
beta4_pH[12,2]	0.581	0.114	0.373	0.573	0.829
beta4_pH[13,2]	0.945	0.045	0.859	0.944	1.037
beta4_pH[14,2]	0.931	0.064	0.803	0.931	1.057
beta4_pH[15,2]	1.011	0.025	0.962	1.011	1.061
beta4_pH[16,2]	1.407	0.077	1.259	1.405	1.559
beta4_pH[1,3]	0.288	0.086	0.135	0.283	0.464
beta4_pH[2,3]	0.052	0.017	0.020	0.051	0.087
beta4_pH[3,3]	0.053	0.010	0.034	0.052	0.072
beta4_pH[4,3]	0.173	0.031	0.113	0.172	0.234
beta4_pH[5,3]	0.777	0.954	0.122	0.563	2.769
beta4_pH[6,3]	0.854	0.679	0.199	0.702	2.493
beta4_pH[7,3]	0.879	0.398	0.346	0.806	1.867
beta4_pH[8,3]	0.418	0.194	0.110	0.396	0.852
beta4_pH[9,3]	0.939	0.937	0.177	0.708	2.954
beta4_pH[10,3]	0.061	0.041	0.003	0.056	0.155
beta4_pH[11,3]	0.364	0.016	0.333	0.363	0.396
beta4_pH[12,3]	0.440	0.041	0.361	0.438	0.524
beta4_pH[13,3]	0.839	0.037	0.771	0.838	0.914
beta4_pH[14,3]	0.825	0.066	0.703	0.821	0.959
beta4_pH[15,3]	0.636	0.017	0.604	0.635	0.668
beta4_pH[16,3]	0.638	0.035	0.573	0.637	0.709
beta4_pH[1,4]	0.624	0.330	0.212	0.559	1.445
beta4_pH[2,4]	0.808	0.875	0.123	0.585	2.852
beta4_pH[3,4]	0.850	0.862	0.113	0.607	3.084
beta4_pH[4,4]	0.725	0.449	0.201	0.617	1.829
beta4_pH[5,4]	0.762	0.881	0.119	0.559	2.531
beta4_pH[6,4]	0.941	1.037	0.204	0.745	2.859
beta4_pH[7,4]	0.896	0.439	0.349	0.825	1.897
beta4_pH[8,4]	0.750	0.632	0.165	0.590	2.315
beta4_pH[9,4]	0.763	0.750	0.135	0.569	2.544
beta4_pH[10,4]	1.454	2.910	0.008	0.404	10.628
beta4_pH[11,4]	0.444	0.017	0.412	0.444	0.477
beta4_pH[12,4]	0.562	0.053	0.458	0.561	0.667
beta4_pH[13,4]	0.925	0.036	0.856	0.924	0.996
beta4_pH[14,4]	0.919	0.073	0.778	0.918	1.070
beta4_pH[15,4]	0.713	0.020	0.674	0.713	0.754
beta4_pH[16,4]	0.686	0.042	0.608	0.684	0.773
beta4_pH[1,5]	0.630	0.333	0.214	0.560	1.383
beta4_pH[2,5]	0.813	0.890	0.108	0.592	2.822
beta4_pH[3,5]	0.855	0.978	0.117	0.588	3.166
beta4_pH[4,5]	0.735	0.490	0.209	0.629	1.967
beta4_pH[5,5]	0.761	0.784	0.105	0.559	2.747
beta4_pH[6,5]	0.919	0.991	0.214	0.727	2.653
beta4_pH[7,5]	0.896	0.439	0.340	0.812	1.900
beta4_pH[8,5]	0.746	0.650	0.164	0.603	2.140
beta4_pH[9,5]	0.760	0.765	0.126	0.560	2.687
beta4_pH[10,5]	1.446	2.986	0.007	0.360	11.347
beta4_pH[11,5]	0.602	0.053	0.507	0.600	0.711
beta4_pH[12,5]	0.583	0.164	0.307	0.564	0.951
beta4_pH[13,5]	0.773	0.063	0.661	0.769	0.906
beta4_pH[14,5]	0.842	0.118	0.633	0.835	1.095
beta4_pH[15,5]	1.145	0.072	1.012	1.144	1.290
beta4_pH[16,5]	18.618	4.019	10.580	18.909	24.686
beta0_pelagic[1]	2.269	0.037	2.201	2.268	2.343
beta0_pelagic[2]	1.537	0.007	1.524	1.537	1.551
beta0_pelagic[3]	0.177	0.275	-0.543	0.249	0.519
beta0_pelagic[4]	0.231	0.197	-0.172	0.243	0.574
beta0_pelagic[5]	1.446	0.156	1.156	1.441	1.777
beta0_pelagic[6]	1.354	0.061	1.239	1.353	1.483
beta0_pelagic[7]	1.480	0.039	1.408	1.478	1.558
beta0_pelagic[8]	1.673	0.043	1.591	1.672	1.761
beta0_pelagic[9]	2.794	0.066	2.670	2.793	2.934
beta0_pelagic[10]	2.603	0.031	2.544	2.603	2.664
beta0_pelagic[11]	-0.164	0.181	-0.573	-0.153	0.142
beta0_pelagic[12]	1.801	0.036	1.731	1.800	1.873
beta0_pelagic[13]	-1.008	0.475	-1.713	-1.103	-0.004
beta0_pelagic[14]	-0.273	0.183	-0.714	-0.255	0.033
beta0_pelagic[15]	-0.210	0.034	-0.301	-0.196	-0.181
beta0_pelagic[16]	-0.719	0.455	-1.524	-0.733	0.026
beta1_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[3]	0.785	0.314	0.375	0.716	1.708
beta1_pelagic[4]	0.814	0.218	0.429	0.811	1.263
beta1_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[11]	1.185	0.263	0.711	1.158	1.824
beta1_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[13]	2.810	0.614	1.477	2.891	3.843
beta1_pelagic[14]	1.813	0.474	1.100	1.713	2.969
beta1_pelagic[15]	0.090	0.162	0.000	0.013	0.603
beta1_pelagic[16]	1.858	0.629	0.827	1.878	3.034
beta2_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[3]	6.006	5.950	0.171	4.272	21.267
beta2_pelagic[4]	7.477	6.054	0.408	6.083	22.113
beta2_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[11]	1.661	1.488	0.220	1.166	5.959
beta2_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[13]	0.188	0.092	0.098	0.162	0.443
beta2_pelagic[14]	0.813	0.828	0.199	0.547	3.195
beta2_pelagic[15]	2.639	1.902	0.138	2.257	7.300
beta2_pelagic[16]	0.560	0.875	0.108	0.271	3.745
beta3_pelagic[1]	28.687	3.868	22.287	28.689	35.101
beta3_pelagic[2]	28.753	3.890	22.338	28.572	35.155
beta3_pelagic[3]	27.911	2.682	23.055	28.252	32.184
beta3_pelagic[4]	26.074	1.241	23.266	26.284	28.643
beta3_pelagic[5]	29.312	6.688	18.659	28.966	40.365
beta3_pelagic[6]	29.541	6.650	18.545	29.556	40.475
beta3_pelagic[7]	29.533	6.672	18.537	29.438	40.461
beta3_pelagic[8]	29.454	6.593	18.557	29.231	40.364
beta3_pelagic[9]	29.598	6.659	18.584	29.574	40.441
beta3_pelagic[10]	29.386	6.480	18.580	29.337	40.245
beta3_pelagic[11]	30.042	1.325	27.596	29.974	32.875
beta3_pelagic[12]	30.000	0.000	30.000	30.000	30.000
beta3_pelagic[13]	31.164	3.996	24.723	30.680	38.969
beta3_pelagic[14]	34.470	1.743	31.516	34.239	38.336
beta3_pelagic[15]	33.391	5.906	23.536	33.778	41.723
beta3_pelagic[16]	29.877	2.591	24.912	30.135	34.665
beta0_yellow[1]	-0.483	0.080	-0.638	-0.483	-0.324
beta0_yellow[2]	0.544	0.079	0.368	0.550	0.682
beta0_yellow[3]	-0.220	0.133	-0.495	-0.221	0.013
beta0_yellow[4]	0.877	0.252	0.236	0.939	1.185
beta0_yellow[5]	0.248	0.335	-0.458	0.272	0.834
beta0_yellow[6]	0.822	0.107	0.614	0.825	1.025
beta0_yellow[7]	0.993	0.095	0.809	0.993	1.176
beta0_yellow[8]	1.023	0.092	0.846	1.022	1.205
beta0_yellow[9]	0.460	0.100	0.262	0.459	0.657
beta0_yellow[10]	0.540	0.072	0.399	0.538	0.683
beta0_yellow[11]	-0.194	0.122	-0.460	-0.184	0.020
beta0_yellow[12]	-1.081	0.094	-1.267	-1.079	-0.904
beta0_yellow[13]	-1.362	0.185	-1.788	-1.345	-1.064
beta0_yellow[14]	-0.248	0.142	-0.546	-0.234	-0.011
beta0_yellow[15]	-1.254	0.139	-1.604	-1.236	-1.037
beta0_yellow[16]	-0.244	0.093	-0.428	-0.244	-0.066
beta1_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[2]	1.011	0.231	0.584	1.003	1.474
beta1_yellow[3]	0.483	0.246	0.020	0.497	0.936
beta1_yellow[4]	1.344	0.619	0.637	1.182	3.268
beta1_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[11]	0.270	0.340	0.001	0.089	1.107
beta1_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[13]	1.740	1.517	0.068	1.102	5.703
beta1_yellow[14]	0.737	1.046	0.012	0.439	4.319
beta1_yellow[15]	0.186	0.245	0.011	0.082	0.923
beta1_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[2]	-7.622	5.816	-21.956	-6.300	-0.542
beta2_yellow[3]	-7.581	5.823	-22.272	-6.216	-0.403
beta2_yellow[4]	-5.598	5.918	-20.326	-3.620	-0.124
beta2_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[11]	-6.842	5.237	-20.016	-5.681	-0.411
beta2_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[13]	-5.189	5.301	-18.787	-3.031	-0.277
beta2_yellow[14]	-6.598	4.883	-18.747	-5.884	-0.436
beta2_yellow[15]	-7.500	5.586	-20.823	-6.429	-0.252
beta2_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta3_yellow[1]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[2]	28.853	1.090	27.305	28.661	31.696
beta3_yellow[3]	32.396	3.206	23.371	32.599	38.267
beta3_yellow[4]	30.104	3.289	23.940	29.394	35.691
beta3_yellow[5]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[6]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[7]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[8]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[9]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[10]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[11]	35.198	3.741	30.306	34.173	43.347
beta3_yellow[12]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[13]	32.476	2.545	29.112	32.276	37.578
beta3_yellow[14]	33.832	3.432	29.187	33.497	42.516
beta3_yellow[15]	36.409	3.941	30.243	36.111	43.576
beta3_yellow[16]	35.000	0.000	35.000	35.000	35.000
beta0_black[1]	0.011	0.075	-0.137	0.011	0.159
beta0_black[2]	1.855	0.032	1.791	1.855	1.917
beta0_black[3]	1.239	0.037	1.165	1.240	1.311
beta0_black[4]	2.423	0.048	2.330	2.422	2.520
beta0_black[5]	5.137	2.370	1.875	4.675	11.386
beta0_black[6]	5.259	2.237	2.570	4.700	11.359
beta0_black[7]	4.102	2.094	1.651	3.529	9.949
beta0_black[8]	0.899	0.226	0.484	0.888	1.380
beta0_black[9]	2.808	0.262	2.314	2.803	3.318
beta0_black[10]	1.286	0.050	1.189	1.285	1.386
beta0_black[11]	2.978	0.052	2.876	2.978	3.079
beta0_black[12]	4.913	0.192	4.539	4.912	5.297
beta0_black[13]	0.458	0.337	0.017	0.390	1.228
beta0_black[14]	2.850	0.108	2.642	2.850	3.064
beta0_black[15]	0.721	0.052	0.622	0.721	0.822
beta0_black[16]	3.772	0.073	3.628	3.773	3.915
beta2_black[1]	0.000	0.000	0.000	0.000	0.000
beta2_black[2]	0.000	0.000	0.000	0.000	0.000
beta2_black[3]	0.000	0.000	0.000	0.000	0.000
beta2_black[4]	0.000	0.000	0.000	0.000	0.000
beta2_black[5]	0.000	0.000	0.000	0.000	0.000
beta2_black[6]	0.000	0.000	0.000	0.000	0.000
beta2_black[7]	0.000	0.000	0.000	0.000	0.000
beta2_black[8]	0.000	0.000	0.000	0.000	0.000
beta2_black[9]	0.000	0.000	0.000	0.000	0.000
beta2_black[10]	0.000	0.000	0.000	0.000	0.000
beta2_black[11]	0.000	0.000	0.000	0.000	0.000
beta2_black[12]	0.000	0.000	0.000	0.000	0.000
beta2_black[13]	-2.621	1.809	-6.963	-2.202	-0.266
beta2_black[14]	0.000	0.000	0.000	0.000	0.000
beta2_black[15]	0.000	0.000	0.000	0.000	0.000
beta2_black[16]	0.000	0.000	0.000	0.000	0.000
beta3_black[1]	25.000	0.000	25.000	25.000	25.000
beta3_black[2]	25.000	0.000	25.000	25.000	25.000
beta3_black[3]	25.000	0.000	25.000	25.000	25.000
beta3_black[4]	25.000	0.000	25.000	25.000	25.000
beta3_black[5]	25.000	0.000	25.000	25.000	25.000
beta3_black[6]	25.000	0.000	25.000	25.000	25.000
beta3_black[7]	25.000	0.000	25.000	25.000	25.000
beta3_black[8]	25.000	0.000	25.000	25.000	25.000
beta3_black[9]	25.000	0.000	25.000	25.000	25.000
beta3_black[10]	25.000	0.000	25.000	25.000	25.000
beta3_black[11]	25.000	0.000	25.000	25.000	25.000
beta3_black[12]	25.000	0.000	25.000	25.000	25.000
beta3_black[13]	37.337	3.315	27.221	38.112	40.912
beta3_black[14]	25.000	0.000	25.000	25.000	25.000
beta3_black[15]	25.000	0.000	25.000	25.000	25.000
beta3_black[16]	25.000	0.000	25.000	25.000	25.000
beta4_black[1]	-0.328	0.092	-0.508	-0.327	-0.146
beta4_black[2]	0.257	0.053	0.155	0.257	0.363
beta4_black[3]	-0.720	0.064	-0.842	-0.721	-0.598
beta4_black[4]	0.636	0.129	0.388	0.634	0.898
beta4_black[5]	0.759	1.585	-1.141	0.455	4.313
beta4_black[6]	0.850	1.474	-0.884	0.533	4.527
beta4_black[7]	0.768	1.460	-0.951	0.490	3.983
beta4_black[8]	0.078	0.297	-0.564	0.091	0.655
beta4_black[9]	1.031	0.821	-0.054	0.858	3.183
beta4_black[10]	0.076	0.066	-0.057	0.077	0.203
beta4_black[11]	-0.609	0.066	-0.738	-0.609	-0.481
beta4_black[12]	1.243	0.300	0.679	1.237	1.856
beta4_black[13]	-1.005	0.077	-1.164	-1.004	-0.858
beta4_black[14]	0.222	0.184	-0.138	0.221	0.592
beta4_black[15]	0.003	0.047	-0.089	0.003	0.093
beta4_black[16]	-0.416	0.113	-0.634	-0.417	-0.192
beta0_dsr[11]	-0.329	0.009	-0.347	-0.329	-0.311
beta0_dsr[12]	1.408	0.043	1.326	1.406	1.495
beta0_dsr[13]	0.155	0.021	0.113	0.155	0.198
beta0_dsr[14]	-0.127	0.017	-0.159	-0.127	-0.094
beta0_dsr[15]	0.205	0.014	0.179	0.205	0.232
beta0_dsr[16]	0.157	0.016	0.127	0.157	0.188
beta1_dsr[11]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[12]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[13]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[14]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[15]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[16]	0.000	0.000	0.000	0.000	0.000
beta4_dsr[11]	0.543	0.035	0.473	0.543	0.613
beta4_dsr[12]	0.808	0.226	0.385	0.803	1.268
beta4_dsr[13]	-0.078	0.041	-0.159	-0.077	0.000
beta4_dsr[14]	-0.208	0.071	-0.348	-0.208	-0.067
beta4_dsr[15]	0.217	0.022	0.175	0.218	0.260
beta4_dsr[16]	-0.047	0.045	-0.136	-0.047	0.041
beta0_slope[11]	-2.209	0.021	-2.251	-2.208	-2.168
beta0_slope[12]	-4.985	0.188	-5.379	-4.980	-4.626
beta0_slope[13]	-1.501	0.034	-1.567	-1.501	-1.437
beta0_slope[14]	-2.712	0.053	-2.822	-2.710	-2.613
beta0_slope[15]	-1.592	0.022	-1.635	-1.592	-1.551
beta0_slope[16]	-2.653	0.042	-2.734	-2.655	-2.569
beta1_slope[11]	0.000	0.000	0.000	0.000	0.000
beta1_slope[12]	0.000	0.000	0.000	0.000	0.000
beta1_slope[13]	0.000	0.000	0.000	0.000	0.000
beta1_slope[14]	0.000	0.000	0.000	0.000	0.000
beta1_slope[15]	0.000	0.000	0.000	0.000	0.000
beta1_slope[16]	0.000	0.000	0.000	0.000	0.000
beta4_slope[11]	-0.335	0.075	-0.490	-0.335	-0.191
beta4_slope[12]	-4.667	0.793	-6.327	-4.606	-3.266
beta4_slope[13]	0.328	0.054	0.225	0.328	0.431
beta4_slope[14]	-0.158	0.152	-0.456	-0.156	0.137
beta4_slope[15]	-0.360	0.033	-0.424	-0.360	-0.294
beta4_slope[16]	-0.139	0.102	-0.344	-0.136	0.053
beta5_slope[11]	3.069	0.107	2.868	3.063	3.291
beta5_slope[12]	2.008	0.369	1.325	1.995	2.762
beta5_slope[13]	0.881	0.079	0.729	0.879	1.038
beta5_slope[14]	3.761	0.365	3.099	3.740	4.543
beta5_slope[15]	1.310	0.073	1.168	1.308	1.456
beta5_slope[16]	2.731	0.230	2.304	2.716	3.206
tau_comp[1]	3.780	0.295	3.221	3.778	4.381
tau_comp[2]	2.859	0.257	2.382	2.852	3.403
tau_comp[3]	0.869	0.074	0.733	0.866	1.020
tau_comp[4]	1.779	0.274	1.304	1.758	2.377
tau_comp[5]	1.262	0.183	0.941	1.249	1.654
sigma_H[1]	0.205	0.048	0.121	0.203	0.308
sigma_H[2]	0.175	0.029	0.127	0.172	0.237
sigma_H[3]	0.188	0.040	0.115	0.186	0.272
sigma_H[4]	0.462	0.077	0.330	0.454	0.636
sigma_H[5]	0.956	0.206	0.582	0.947	1.399
sigma_H[6]	0.510	0.188	0.139	0.501	0.907
sigma_H[7]	0.312	0.062	0.214	0.305	0.454
sigma_H[8]	0.425	0.087	0.291	0.415	0.621
sigma_H[9]	0.516	0.117	0.328	0.502	0.778
sigma_H[10]	0.209	0.039	0.144	0.206	0.294
sigma_H[11]	0.268	0.044	0.193	0.263	0.364
sigma_H[12]	0.381	0.144	0.197	0.337	0.725
sigma_H[13]	0.218	0.037	0.156	0.215	0.299
sigma_H[14]	0.499	0.088	0.346	0.493	0.689
sigma_H[15]	0.231	0.038	0.168	0.226	0.318
sigma_H[16]	0.216	0.040	0.148	0.212	0.306
lambda_H[1]	2.186	3.740	0.080	0.998	11.214
lambda_H[2]	5.024	4.537	0.504	3.743	17.014
lambda_H[3]	5.305	8.178	0.246	2.695	26.661
lambda_H[4]	0.006	0.004	0.001	0.005	0.017
lambda_H[5]	3.367	10.883	0.027	0.650	24.778
lambda_H[6]	3.837	9.178	0.008	0.117	32.088
lambda_H[7]	0.012	0.008	0.002	0.010	0.035
lambda_H[8]	7.521	10.109	0.077	3.933	36.315
lambda_H[9]	0.015	0.010	0.003	0.013	0.041
lambda_H[10]	0.229	0.330	0.025	0.153	0.843
lambda_H[11]	0.237	0.362	0.011	0.106	1.167
lambda_H[12]	3.830	5.103	0.184	2.207	17.023
lambda_H[13]	3.076	2.762	0.254	2.328	10.405
lambda_H[14]	3.320	4.185	0.214	2.048	13.862
lambda_H[15]	0.019	0.047	0.003	0.014	0.057
lambda_H[16]	0.586	1.034	0.035	0.282	3.052
mu_lambda_H[1]	4.052	1.888	1.038	3.844	8.116
mu_lambda_H[2]	3.499	1.928	0.440	3.282	7.811
mu_lambda_H[3]	3.394	1.810	0.739	3.140	7.545
sigma_lambda_H[1]	8.032	4.321	1.743	7.277	17.824
sigma_lambda_H[2]	7.789	4.734	0.760	6.959	18.246
sigma_lambda_H[3]	6.174	3.903	0.982	5.399	15.906
beta_H[1,1]	6.501	1.272	3.378	6.708	8.454
beta_H[2,1]	9.778	0.602	8.497	9.814	10.901
beta_H[3,1]	7.819	0.824	5.779	7.937	9.161
beta_H[4,1]	10.442	7.486	-4.289	10.427	24.719
beta_H[5,1]	-0.319	2.561	-5.577	-0.165	4.194
beta_H[6,1]	2.539	4.490	-7.827	3.836	8.280
beta_H[7,1]	1.142	5.676	-10.980	1.550	11.348
beta_H[8,1]	1.474	4.662	-2.705	1.233	3.865
beta_H[9,1]	13.665	5.655	2.708	13.545	24.920
beta_H[10,1]	6.865	1.866	2.977	6.928	10.557
beta_H[11,1]	4.911	3.537	-2.671	5.620	9.994
beta_H[12,1]	2.596	1.082	0.659	2.524	4.934
beta_H[13,1]	9.105	0.891	7.220	9.167	10.589
beta_H[14,1]	2.141	1.028	0.134	2.146	4.147
beta_H[15,1]	-7.108	3.356	-13.199	-7.305	-0.095
beta_H[16,1]	3.839	2.900	-0.723	3.401	10.311
beta_H[1,2]	7.778	0.268	7.232	7.785	8.283
beta_H[2,2]	9.880	0.165	9.556	9.879	10.197
beta_H[3,2]	8.799	0.207	8.385	8.803	9.216
beta_H[4,2]	3.030	1.419	0.289	3.018	5.817
beta_H[5,2]	1.736	0.991	-0.194	1.736	3.655
beta_H[6,2]	5.854	1.136	3.350	6.009	7.682
beta_H[7,2]	2.418	1.059	0.464	2.364	4.619
beta_H[8,2]	2.956	1.277	1.063	3.162	4.292
beta_H[9,2]	3.171	1.078	1.111	3.137	5.371
beta_H[10,2]	8.025	0.371	7.250	8.035	8.711
beta_H[11,2]	9.821	0.636	8.867	9.706	11.210
beta_H[12,2]	3.941	0.375	3.262	3.924	4.718
beta_H[13,2]	9.157	0.246	8.702	9.148	9.654
beta_H[14,2]	4.028	0.339	3.389	4.024	4.719
beta_H[15,2]	11.566	0.606	10.307	11.579	12.727
beta_H[16,2]	4.492	0.823	2.998	4.462	6.213
beta_H[1,3]	8.462	0.262	7.985	8.455	8.989
beta_H[2,3]	9.854	0.150	9.562	9.853	10.149
beta_H[3,3]	9.459	0.168	9.141	9.453	9.808
beta_H[4,3]	-2.199	0.855	-3.909	-2.194	-0.580
beta_H[5,3]	3.906	0.595	2.679	3.911	5.071
beta_H[6,3]	8.543	1.201	6.710	8.375	10.916
beta_H[7,3]	-2.572	0.634	-3.804	-2.563	-1.336
beta_H[8,3]	5.329	0.587	4.695	5.229	6.610
beta_H[9,3]	-2.696	0.697	-4.121	-2.677	-1.401
beta_H[10,3]	8.560	0.278	8.016	8.564	9.111
beta_H[11,3]	8.502	0.293	7.895	8.521	9.005
beta_H[12,3]	5.272	0.321	4.528	5.313	5.776
beta_H[13,3]	8.849	0.176	8.491	8.854	9.180
beta_H[14,3]	5.785	0.267	5.213	5.807	6.248
beta_H[15,3]	10.277	0.301	9.696	10.274	10.887
beta_H[16,3]	5.986	0.614	4.741	6.001	7.120
beta_H[1,4]	8.096	0.187	7.710	8.109	8.432
beta_H[2,4]	9.844	0.130	9.585	9.850	10.080
beta_H[3,4]	9.940	0.156	9.594	9.956	10.205
beta_H[4,4]	11.815	0.468	10.884	11.809	12.715
beta_H[5,4]	5.596	0.751	4.360	5.519	7.285
beta_H[6,4]	7.027	0.943	5.191	7.114	8.501
beta_H[7,4]	8.547	0.354	7.874	8.545	9.258
beta_H[8,4]	6.707	0.282	6.126	6.731	7.136
beta_H[9,4]	7.423	0.442	6.554	7.420	8.286
beta_H[10,4]	7.519	0.224	7.089	7.512	7.992
beta_H[11,4]	9.504	0.195	9.132	9.505	9.878
beta_H[12,4]	7.254	0.198	6.880	7.251	7.667
beta_H[13,4]	9.063	0.141	8.781	9.068	9.333
beta_H[14,4]	7.838	0.215	7.436	7.834	8.276
beta_H[15,4]	9.482	0.225	9.046	9.482	9.922
beta_H[16,4]	9.403	0.246	8.940	9.397	9.902
beta_H[1,5]	8.899	0.151	8.579	8.904	9.195
beta_H[2,5]	10.674	0.101	10.484	10.671	10.882
beta_H[3,5]	10.765	0.166	10.470	10.751	11.101
beta_H[4,5]	8.331	0.490	7.369	8.320	9.334
beta_H[5,5]	5.174	0.613	3.754	5.252	6.154
beta_H[6,5]	9.034	0.621	8.048	8.947	10.335
beta_H[7,5]	6.576	0.348	5.896	6.574	7.245
beta_H[8,5]	8.248	0.229	7.886	8.229	8.724
beta_H[9,5]	8.206	0.454	7.312	8.201	9.091
beta_H[10,5]	10.046	0.231	9.591	10.047	10.479
beta_H[11,5]	11.590	0.228	11.141	11.594	12.034
beta_H[12,5]	8.614	0.188	8.240	8.611	8.987
beta_H[13,5]	10.140	0.134	9.883	10.135	10.416
beta_H[14,5]	9.252	0.218	8.844	9.244	9.701
beta_H[15,5]	11.310	0.234	10.843	11.312	11.787
beta_H[16,5]	9.839	0.187	9.450	9.842	10.198
beta_H[1,6]	10.118	0.236	9.737	10.086	10.688
beta_H[2,6]	11.475	0.118	11.236	11.475	11.700
beta_H[3,6]	10.612	0.159	10.277	10.620	10.901
beta_H[4,6]	12.917	0.884	11.135	12.955	14.667
beta_H[5,6]	5.684	0.624	4.579	5.659	7.007
beta_H[6,6]	8.707	0.585	7.435	8.738	9.783
beta_H[7,6]	10.552	0.603	9.403	10.532	11.721
beta_H[8,6]	9.510	0.302	8.984	9.533	9.972
beta_H[9,6]	8.261	0.773	6.764	8.258	9.830
beta_H[10,6]	9.136	0.344	8.376	9.162	9.734
beta_H[11,6]	10.688	0.367	9.894	10.712	11.344
beta_H[12,6]	9.385	0.238	8.967	9.372	9.902
beta_H[13,6]	11.086	0.156	10.781	11.082	11.402
beta_H[14,6]	9.923	0.286	9.331	9.931	10.461
beta_H[15,6]	10.553	0.416	9.726	10.555	11.370
beta_H[16,6]	10.600	0.258	10.052	10.603	11.091
beta_H[1,7]	10.162	1.266	6.833	10.437	11.939
beta_H[2,7]	12.106	0.524	10.988	12.139	13.062
beta_H[3,7]	10.283	0.688	8.783	10.328	11.480
beta_H[4,7]	1.398	4.534	-7.517	1.408	10.208
beta_H[5,7]	5.890	2.036	1.955	5.922	9.889
beta_H[6,7]	6.743	3.355	-1.307	7.741	11.084
beta_H[7,7]	7.674	3.047	1.812	7.697	13.510
beta_H[8,7]	10.927	1.154	9.170	10.873	12.805
beta_H[9,7]	5.018	4.006	-3.178	5.022	12.588
beta_H[10,7]	10.079	1.612	7.314	9.915	13.728
beta_H[11,7]	11.370	1.776	8.262	11.207	15.352
beta_H[12,7]	9.795	0.951	7.707	9.884	11.428
beta_H[13,7]	11.703	0.712	10.044	11.770	12.909
beta_H[14,7]	10.632	0.941	8.656	10.671	12.385
beta_H[15,7]	13.961	2.192	9.607	13.974	18.213
beta_H[16,7]	12.781	1.382	10.471	12.655	15.910
beta0_H[1]	8.879	17.454	-27.803	8.891	47.069
beta0_H[2]	10.382	7.727	-5.007	10.394	26.929
beta0_H[3]	9.256	10.151	-12.403	9.436	29.919
beta0_H[4]	0.734	186.498	-361.080	-0.853	371.596
beta0_H[5]	4.105	29.020	-50.014	4.388	56.644
beta0_H[6]	7.992	57.690	-123.755	7.799	133.766
beta0_H[7]	6.341	135.580	-270.503	5.507	279.325
beta0_H[8]	7.678	38.197	-18.926	6.531	35.883
beta0_H[9]	6.216	122.080	-237.789	6.201	251.139
beta0_H[10]	7.885	38.485	-71.999	8.063	84.549
beta0_H[11]	10.622	53.979	-98.936	10.597	123.651
beta0_H[12]	6.666	12.963	-18.108	6.545	31.515
beta0_H[13]	9.704	11.287	-12.093	9.779	30.364
beta0_H[14]	7.331	11.413	-15.018	7.142	30.301
beta0_H[15]	10.111	112.798	-224.869	11.319	227.900
beta0_H[16]	9.284	29.495	-53.377	8.952	72.667